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ABSTRACT 

We extend the formalism of the relativistic thermal and kinematic Sunyaev-Zel'dovich effects 
and include the polarization of the cosmic microwave background photons. We consider the 
situation that a cluster of galaxies is moving with a velocity /3 = u/c with respect to the cosmic 
microwave background radiation. In the present formalism, polarization of the scattered cosmic 
microwave background radiation caused by the proper motion of a cluster of galaxies is naturally 
derived as a special case of the kinematic Sunyaev-Zel'dovich effect. The relativistic corrections 
are included also in a natural way. Our results are in complete agreement with the recent results 
of the relativistic corrections obtained by Challinor, Ford & Lasenby with an entirely different 
method as well as the nonrelativistic limit obtained by Sunyaev & Zel'dovich. The relativistic 
correction becomes to be significant in the Wien region. 

Subject headings: cosmic microwave background cosmology: theory — galaxies: clusters: 
general — radiation mechanisms: thermal — polarization — relativity 
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1. INTRODUCTION 

The present authors (Itoh, Kohyama, & Nozawa 1998; Nozawa, Itoh, & Kohyama 1998) have recently 
given accurate relativistic corrections to the thermal and kinematic Sunyaev-Zel'dovich effects (Zel'dovich 
& Sunyaev 1969; Sunyaev & Zel'dovich 1972,1980a, 1980b, 1981). Their method was based on the kinetic 
equation for the photon distribution function using a relativistically covariant formahsm (Berestetskii, 
Lifshitz, & Pitaevskii 1982; Buchler & Yueh 1976). By using the generaUzed Kompaneets equation 
(Kompaneets 1957; Weymann 1965), they have derived analytic expressions for the thermal and kinematic 
Sunyaev-Zel'dovich effects as a power series of 9^ = kBTe/mc? and /?, where Tg is the electron temperature 
and /? is the peculiar velocity of the cluster divided by the velocity of light. It has been shown by several 
groups (Stebbins 1997; Challinor & Lasenby 1998; Itoh, Kohyama, & Nozawa 1998) that the results for the 
thermal Sunyaev-Zel'dovich effect obtained by the power series expansion agree with the previous numerical 
calculations by RephaeH (1995) and by RephaeU & Yankovitch (1997), thereby proving the validity of their 
method. In particular, the convergence of the power scries expansion has been carefully studied in the paper 
by Itoh, Kohyama, & Nozawa (1998), whore the analytic expressions up to 0{6"l) have been derived and 
the results have been compared with those of the direct numerical integration of the Boltzmann equation. 
It has been shown that the power series expansion approximation is sufficiently accurate for the region 
ksTe < 15keV by taking into account up to 0{9l) contributions. 

In the paper by Nozawa, Itoh, & Kohyama (1998), a formalism for the kinematic Sunyaev-Zel'dovich 
effect for the clusters of galaxies with a peculiar velocity /? has been derived by applying the Lorentz 
boost to the standard formalism of the extended Kompaneets equation. With the power series expansion 
approximation in terms of the electron temperature 6e and the peculiar velocity /3, an analytic expression 
for the kinematic Sunyaev-Zel'dovich effect which includes the relativistic corrections of up to 0{(i9l) and 
0{fP'9e) has been derived. It has been found that the relativistic correction is significant. Similar works 
have been reported by Sazonov & Sunyaev (1998), and also by Challinor & Lasenby (1999). They are in 
essential agreement with the works of the present authors. 

In the present paper we address ourselves to the problem of the polarization of the cosmic microwave 
background radiation (CMBR) caused by the proper motion of the cluster of galaxies. This problem has been 
studied by Sunyaev & Zel'dovich (1980b) in the nonrelativistic limit. The polarization Sunyaev-Zel'dovich 
effect enables the measurement of the tangential velocities of the clusters of galaxies, thus contributing to 
the measurement of the radial velocities complemented by the measurement of the velocities along the line 
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of sight which is achieved by means of the observation of the kinematic Sunyaev-Zel'dovich effect. There 
exist, however, competing polarization effects connected with the finite depth of the intracluster plasma. 
These effects have been discussed by Sunyaev & Zel'dovich (1980b), and more recently by Sazonov & 
Sunyaev (1999). In the present paper we solve this problem by extending our already-established formalism 
to the case of polarized photons. With the present formalism the effect of the polarization of the CMBR 
can be derived on the same footing as the calculation of the thermal and kinematic Sunyaev-Zel'dovich 
effects. We thereby take into account the relativistic corrections accurately. 

After the submittal of the original manuscript of the present paper to The Astrophysical Journal, we 
became aware of the similar works by Audit & Simmons (1998) as well as by Hansen & Lilje (1999). More 
recently the papers by Sazonov & Sunyaev (1999) and by Challinor, Ford & Lasenby (1999) have been 
submitted on the similar subject (the polarization of the CMBR). Through the communication with these 
authors, we have benefited a lot and have revised the present paper. In particular. Dr. Challinor has kindly 
pointed out an error in the original manuscript of the present paper. In the present revised manuscript, we 
have fully taken into account their comments and have corrected the error in the original manuscript. With 
the present formalism we have confirmed the recent result obtained by Challinor, Ford & Lasenby (1999) 
with an entirely different method as well as the nonrelativistic limit obtained by Sunyaev & Zel'dovich 
(1980b). 

The present paper is organized as follows: the Kompaneets equation for polarized photons will be 
derived in § 2. Results of the calculation will be presented in § 3. Concluding remarks will be given in § 4. 

2. LORENTZ BOOSTED KOMPANEETS EQUATION 

In the present section we will extend the Kompaneets equation studied by Nozawa, Itoh, & Kohyama 
(1998) to the case of polarized photons. We will consider a cluster of galaxies moving with a peculiar 
velocity with respect to the cosmic microwave background radiation (CMBR). We wiU formulate the kinetic 
equation for the photon distribution function using a relativistically covariant formalism (Berestetskii, 

Lifshitz, & Pitacvskii 1982; Buchlcr & Yuch 1976). As a reference system, we choose the system which is 
fixed to CMBR. The z-axis is fixed to the line connecting the cluster of galaxies (CG) and the observer. 
(We assume that the observer is fixed to the CMBR frame.) We fix the positive direction of the z-axis as 
the direction of the propagation of a photon from the cluster to the observer. In this reference system, the 
center of mass of CG is moving with a peculiar velocity /?(= v/c) with respect to the CMBR. For simplicity. 
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we choose the direction of the velocity in the x-z plane, i.e. (3 = {(3x, 0, pz)- 



In order to describe the Compton scattering of polarized photons by unpolarized electrons, we derive 
the Stokes parameters using the polarization density matrix (Berestetskii, Lifshitz, & Pitaevskii (1982)). 
We here emphasize the importance of using the relativistically covariant formalism. Sunyaev & Zel'dovich 
(1980b) have shown that the degree of the polarization of the cosmic microwave background photon caused 
by the proper motion of the cluster of galaxies is of order On the other hand, Nozawa, Itoh, & Kohyama 
(1998) have derived the relativistic corrections to the kinematical Sunyaev-Zel'dovich effect which are of 
order /3^. Therefore, one has to be extremely careful to take into account all the relevant terms of order 
when one calculates the polarization Sunyaev-Zel'dovich effect. This is the reason why we adopt the 
relativistically covariant formalism in this paper. 

In the CMBR frame, the time evolution of the photon distribution n{Lo) is written as follows. 



X { Xk+p^k'+p' n{w)[l + n{u')]f{E) - Xw+p'^k+p n{w')[l + n{w)]f{E') } , 



(2.1) 



where k, k' are the four- momenta of photons and p,p' are the four-momenta of electrons, respectively. In 
equation (2.1) the first term corresponds to the Compton scattering k + p ^ k' + p'. The explicit form is 
given as follows. 



^k+p—^k'+p' 



X 1 + fj-c 



Xi) 



(2.2) 



where 



X = 
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X = (s — m^) I'i 
y = [m^ - u) / 



(2.7) 
(2.8) 
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/\2 



s = {p + ky = (p' + k') 

u = {p- k'f = ip' - kf 



(2.9) 
(2.10) 



In equation (2.2), rf = (77i,r?2,%) is the Stokes parameter which corresponds to the selection of the 
polarization of the photon of momentum k by the detector (Berestetskii, Lifshitz, & Pitaevskii (1982)). In 
equations (2.4)-(2.6), e^^^ and e^^' are the polarization vectors for photons of the momentum k. We choose 
the explicit forms as follows: 



k = (0,0,w), 
kxk' 



k = {lo, k) , 
eW = (0,e«), e« 

I ft X r 

e(2) = (0,e(2)), g(2) ^ ft X (ft X ftp 



ft X (ft X ft') 



^k'+p'—i-k+p — X 



where 
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1 8 
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= (1) .«'p(2) 



gv-y . p- gv-y .p _ I 1 _ g(l) . p' g(2) . 

X yj \x yj X I 



1 1\ 1 



1 1 
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'"■"f -( 



e(2).p' 



X y J X x^V /J 



Inserting equations (2.2) and (2.14) into equation (2.1), we obtain 

dn{uj) _ ^ [ <fp ^3/^3,/ 



dt 



-2 



2 j ^d^p'd^k'W {n{u;)[l + niio'mE) - niio')[l + n{ujmE')} 
"^'^ d^p'd^k'Wrf-\^ C^^n{uj)[l + niuj')]fiE) 



(27r)3 



where 



^ (eV47r)^<54(p^fe_p/_fe,) 
2a;a;'i;£' 



(2.11) 
(2.12) 

(2.13) 



Similarly, the second term in equation (2.1) corresponds to the Compton scattering ft' +p' — > ft + p. 
The explicit form is as follows: 



(2.14) 

(2.15) 
(2.16) 

(2.17) 



(2.18) 



(2.19) 
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is the transition probability corresponding to the Compton scattering of an unpolarized photon by an 
unpolarized electron. In equation (2.18) the first line corresponds to the unpolarized case, which has been 
studied in Itoh, Kohyama, & Nozawa (1998) and in Nozawa, Itoh, & Kohyama (1998). The second and 
third lines correspond to the polarized photon contribution. 



Equation (2.18) is further simplified. With the help of the following relations 



it is straightforward to show that 



e(2) . (p + fc - fc') = e(2)-p - e(2).fc', 



(2.20) 

(2.21) 

(2.22) 



Therefore we have 



dn{w) 
dt 



■■I 



{2n) 
(27r)3 



d^p'd^k'Wr]-C{n{co)[l + n{u)')]f{E) 

-n{io')[l+n{u,)]f{E')} . 



(2.23) 



The vector ( represents the Stokes parameter of the photon k which is polarized due to the Compton 
scattering (Berestetskii, Lifshitz, & Pitaevskii (1982)). 

Next we transform the Stokes parameter C=(CijC2,C3) to the Stokes parameter ^=(^,^2,^3) defined 
with respect to the coordinate system xyz fixed to the CMBR. According to Landau & Lifshitz (1975), we 
obtain 



6 = Ci cos 2 
6 = C2 , 
6 = Ci sin 2 



Ca sin 2 



(■3 cos 2 



(2.24) 

(2.25) 
(2.26) 



where (f)'^ is the azimuthal angle of the vector k' with respect to the x-z plane. Since C2 = 0, we obtain 
^2 = 0. Therefore, we have 

- cos 2(j)'^ - sin 2(j)'^ 
sin 2(/)!j, — cos 2(/>^y 

where = and Cs = are given by equations (2.15) and (2.17) 





(2.27) 
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Thus, the time evolution of the Stokes parameter corresponding to the photon polarization in the 
CMBR frame is written as follows (Acquista & Anderson 1974; Nagirner 1994): 



d_ 
di 



n{uj) 



^3 



-d'^p'dPk'W 



-cos 20^ -s\n2(t)'^ 
sva-2(j)'^ -cos 20!^ 



Ci 

C3 



X { n{Lo)[l + n{J)]f{E) - n{^')[l + n{uj)]f{E') } . 



(2.28) 



The electron Fermi distribution functions in the initial and final states are defined in the CG frame. They 
are related to the electron Fermi distribution functions in the CMBR frame as follows (Landau & Lifshitz 
1975): 



f{E) = MEc), 

f{E') = fciE'c), 

Ec = j[e-P-p^ , 

E'c = j^E'-P-if) 
1 

7 = 



(2.29) 
(2.30) 
(2.31) 
(2.32) 
(2.33) 



where the sufBx C denotes the CG frame. We caution the reader that the electron distribution function 
is anisotropic in the CMBR frame for P 0. By ignoring the degeneracy effects, we have the relativistic 
Maxwellian distribution for electrons with temperature Te as follows: 

-1 



fc{Ec) 



g{(Ec-m)-(nc-m)}/kBTe _j_ 
-{{Ec-m)-{nc-m)}/kBn 



(2.34) 



where {fic — fn) is the non-relativistic chemical potential of the electron measured in the CG frame. We 
now introduce the quantities 



Ax = 



ksTe ' 
iv' — Ul 



(2.35) 
(2.36) 



Substituting equations (2.29) - (2.36) into equation (2.28), we obtain 



d_ 
dt 



n{uj) 



6 



= -2 



Td^p'd^k'W 



(27r)3 

X fc{Ec)[{l+n{uj')}n{io) 



cos 20^ 
sin 2^; 



-sin 20; 
-cos 2(j)'^ 



Ci 

C3 



{l + n{uj)}n{ijj') e 



(2.37) 



-9- 



where k and k' are the unit vectors in the directions of k and k', respectively. Equation (2.37) is our basic 
equation. 

3. RESULTS OF THE CALCULATION 



We expand the right-hand-side of equation (2.37) in powers of Aa; by assuming Ax <C 1 as has been 
done by Nozawa, Itoh, & Kohyama (1998): 



d_ 
dt 



n(uj) 



where 



6 
6 



— /i,o + n{l + n)Ii^i 



, dn 



+ 2 
+ • • • 

+ 2n 



hfl + 2(1 + n)— h,i + n{l + n) 
^ 73,0 + 3(1 + n)-^h,i + 3(1 + n)—l3,2 + n(l + n) 73,3 



(1 + „)J„ + _J, + _J, + _J3, 



-cos 20' -sin 20' 



^'■J W \^ gi^20; -cos 20; 

X fc{Ec) {Axf e=^-<^<k-k')^^ (1 - /3 • fc')^ , 



Ci 

C3 



(3.1) 



(3.2) 



Jk 



k\j 



^■J (27r)3 \^ gi^20; -cos 20; y V Cs 



xfc{Ec){Axf[l - e^^'^^-f^-^')) 



(3.3) 



For most of the clusters of galaxies, /3 <C 1 is realized. For example, (3 » 1/300 for a typical value of the 
peculiar velocity v = l,000km/s. Therefore it should be sufficient to expand in powers of (3 and retain up 
to O{0^) contributions. We assume the initial photon distribution of the CMBR to be Planckian with a 
temperature Tq: 



no(X) = 



where 



X = 



- 1 



ksTo 



(3.4) 



(3.5) 
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Assuming Tq/T,, <C 1, we obtain 
1 



no{X) 



A 



n{X) 




fix 

y 



yXe 



X 



'^^(3l{Fo + e,F,+elF^) 



P sin O-y , 
ar j diNe , 



(3.6) 

(3.7) 
(3.8) 

(3.9) 



° 20 ' 



(3.10) 



3 21 



:x- — [x^ + -S^]+ — (x^ + 2XS' 

4 5V 2 y 280 V 



867 



(3.12) 



X 

s 



Xcoth f — 



sinh ( — 
2 



(3.13) 
(3.14) 



where 6j is the angle between the directions of the pecuHar velocity (/?) and the photon momentum (fc) 
which is chosen as the 2:-direction. A^e is the electron number density in the CG frame. The integral in 
equation (3.8) is over the photon path length in the cluster, and ctt is the Thomson cross section. 

Equations (3.6) (3.14) are in complete agreement with the recent result obtained by Challinor, Ford 
& Lasenby (1999) with an entirely different method. It should be emphasized that the effect of the 
polarization of the CMBR has been derived on the same footing as the calculation of the thermal and 
kinematic Sunyaev-Zel'dovich effects by using the present formalism. It is also worthwhile to mention here 
that the polarization (equation (3.6)) is proportional to 0^, where (3x is the transverse component of the 
peculiar velocity of CG. In the present case the polarization of the CMBR is nonzero, because the electron 
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distribution is anisotropic in the CMBR frame for the nonzero pecuhar velocity. It should be also noted 
that we have neglected higher order relativistic corrections such as 0{6l) terms in deriving equation (3.6). 
This is because the observation of the higher order term in the polarization of the CMBR caused by the 
motion of the CG will not be feasible in the near future. Therefore it should be sufficient to neglect higher 
order relativistic corrections. 

In the Rayleigh-Jeans limit X ^ with 9e = 0, equation (3.6) gives the polarization degree 

P = ^yf^l, (3-15) 

which agrees with the result of Sunyaev & Zel'dovich (1980b). Similar results have been reported in the 
more recent works: Audit & Simmons (1998); Hansen & Lilje (1999); Sazonov & Sunyaev (1999); Challinor, 
Ford & Lasenby (1999). 

In Figure 1 we show the graph of the function 

^ _ 1 A[n(X)6] _ Xe^ ( p ^ f) F ^ ft^F \ C^^f^^ 

It is clear from Figure 1 that the value of the polarization function E is positive (linear polarization in 
the plane formed by the velocity vector of the cluster and the vector which connects the cluster and the 
observer) and small for the Rayleigh-Jeans region, however, it quickly becomes large for the Wien region, 
where the coefficient is larger than 10 at X > 11. Finally we define the distribution of the spectral intensity 
of the polarized radiation as follows: 

AIpoi = X'A[n{X)C3] = ypl-^p^ ( Fo + e,F^ + BIf^ ) . (3.17) 

(e^ - 1) 

The graph of the function Alpd/y is shown in Figure 2. It has been found that the function has a maximum 
value of 11 X 10~^ at X w 5 for UbT = lOkeV, /3x = 1/300. It is clear that the relativistic correction 
becomes significant for large values of X, typically for X > 5. 



4. CONCLUDING REMARKS 

We have calculated the Stokes parameter corresponding to the polarization of the CMBR caused 
by the proper motion of a cluster of galaxies. The calculation is based on the covariant formalism. 
With the present formalism the polarization of the CMBR is calculated as a natural extension of the 
kinematic Sunyaev- Zel'dovich effect calculated by Nozawa, Itoh, & Kohyama (1998). Thus the relativistic 
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corrections have been derived from first principles. We have confirmed the recent result of the relativistic 
corrections obtained by Challinor, Ford, & Lasenby (1999) with an entirely different method as well as 
the nonrelativistic limit obtained by Sunyaev & Zel'dovich (1980b). The relativistic correction becomes 
significant for large values of X, X > 5. The distribution of the spectral intensity of the polarized radiation 
has been also calculated. The maximum value of Alpoi/y is 11 x 10~^ at X w 5 for ksT = lOkeV, 
f3x = 1/300. At present the observation of the polarization of the CMBR caused by the motion of a cluster 
of galaxies is not feasible. However, it is hoped that its observation becomes possible in the future. (See 
Sazonov & Sunyaev (1999) for the future projects.) 
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error in the original manuscript of the present paper. We also wish to thank the referee for many valuable 
suggestions on the original manuscript which have helped us tremendously in the revision of the paper. This 
work is financially supported in part by the Grant-in- Aid of the Japanese Ministry of Education, Science, 
Sports, and Culture under the contract #08304026 and #10640289. 
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Figure Captions 



• Fig.l. Graph of the function S defined by equation (3.16) for ksTg = lOkeV. The dotted curve is the 
leading order (Fq term) contribution. The dashed curve is the contribution which includes the 0{9e) 
term. The solid curve is the full contribution which includes the 0{6l) term. 

• Fig.2. Graph of Mpoi/y defined by equation (3.17) for the case of fcsTe = lOkeV, = 1/300. 
The dashed curve is the contribution which includes the 0{9e) term. The solid curve is the full 
contribution which includes the 0(9'^) term. 



